Hyperfine interactions in two dimensional HgTe topological insulators 
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We study the hyperfine interaction between the nuclear spins and the electrons in a HgTe quantum 
well, which is the prime experimentally realized example of a two-dimensional topological insulator. 
The hyperfine interaction is a naturally present, internal source of broken time-reversal symmetry 
from the point of view of the electrons. The HgTe quantum well is described by the so-called 
Bernevig-Hughes-Zhang (BHZ) model. The basis states of the BHZ model are combinations of 
both S- and P-like symmetry states, which means that three kinds of hyperfine interactions play 
a role: (i) The Fermi contact interaction, (ii) the dipole-dipole like coupling and (iii) the electron 
orbital to nuclear spin coupling. We provide benchmark results for the forms and magnitudes of 
these hyperfine interactions within the BHZ model, which give a good starting point for evaluating 
hyperfine interactions in any HgTe nanostructure. We apply our results to the helical edge states of a 
HgTe two-dimensional topological insulator and show how their total hyperfine interaction becomes 
anisotropic and dependent on the orientation of the sample edge within the plane. Moreover, for the 
helical edge states the hyperfine interactions due to the P-like states can dominate over the S-like 
contribution in certain circumstances. 



I. INTRODUCTION 



A topological insulator (TI) host gapless surface or 
edge states, while the bulk of the material has an insulat- 
ing energy gap.— — In three-dimensional TIs the gapless 
surface states are spin-polarized two-dimensional (2D) 
Dirac fermions, whereas 2D TIs contain one-dimensional 
(ID) helical edge states. The helical edge states appear 
in counterpropagating pairs, and the states with equal 
energy and opposite wave numbers, k and —k, form a 
Kramers pair. Thus, elastic scattering from one helical 
edge state (HES) to the other one within a pair cannot be 
induced by time-reversal invariant potentials e.g. stem- 
ming from impurities.— Therefore, the transport through 
a 2D TI is to a large extend ballistic with a quantized 
conductance of e 2 /h per pair of HESs. This highlights 
the central role of time-reversal symmetry in TIs. 

Quantized conductance have recently been measured in 
micron-sized samples in HgTe quantum wellsj£~— which 
to date is the most important experimental demonstra- 
tion of a 2D TI. Evidence of edge state transport was 
found in both two-terminal 6 - and multi-terminal^ devices. 
Moreover, clever experiments combining the metallic spin 
Hall effect and a 2D TI in a HgTe quantum well (QW) 
demonstrated the connection between the spin and prop- 
agation direction* 1 ^ However, also deviations from per- 
fect conductance have been observed in longer HgTe 
devices j & 7 ' n i 12 which could stem from e.g. inelastic scat- 
tering mechanisms J^r— The effect of external magnetic 
fields have also been considered! 6 ' 8 ' 19 " — The TI state 
in HgTe QWs was predicted by Bernevig, Hughes and 
Zhang (BHZ)2£ by using a simplified k • p model con- 
taining states with S- and P-like symmetry, respec- 
tively. They found that beyond a critical thickness of 
the HgTe QW, the TI state would appear as confirmed 
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FIG. 1. Hyperfine coupling between the BHZ basis states 
{\H+), \E+), \E-), \H-)}, which have the total angular 
momentum projections on the z-axis mj as indicated. Every 
increase (decrease) of irij in the electronic sector is accompa- 
nied by a decrease (increase) of a nuclear spin due to angular 
momentum conservation. Hyperfine interactions due to both 
S- and P-like states connect the two time-reversed blocks of 
the BHZ model (red full arrows). However, only hyperfine in- 
teractions due to P-like states connect states within a single 
time-reversed block (blue dashed arrows). 



experimentally^— Furthermore, interesting experimen- 
tal progress on 2D TI properties has also be achieved in 
InAs/GaSb QWs^ 6 - as proposed theoretically^ 

Hyperfine (HF) interactions between the electron and 
nuclear spins can play an important role in nanostruc- 
tures - even though it is often weak<22r— For instance in 
quantum dots, HF interactions can limit the coherence 
of single electronic spins^ 2 - - — and, moreover, it can even 
lead to current hysteresis due to bistability of the dynam- 
ical nuclear spin polarization.— ~— HF-induced nuclear 
spin ordering in interacting I D 40 ' 41 and 2 D 42 ' 43 systems 
have also been considered theoretically. Most studies 



consider the so-called contact HF interaction^ - — which 
is relevant for electrons in orbital states with S'-like sym- 
metry e.g. the conduction band in GaAs. However, for 
P-like orbital states - such as the valence band in GaAs - 
the contact HF interaction is absent. Nevertheless, other 
anisotropic HF interactions are present for P-like states 
such as the dipole-dipolc-likc HF interaction^ - — which 
can play a significant role e.g. for the decoherence of a 
hole confined in a quantum doti 44 ' 45 ' 48 

HF interactions and dynamical nuclear spin polariza- 
tion have also been investigated in the context of in- 
teger quantum Hall systems^ 9 - - — which contain unidi- 
rectional edge states. Here HF-induced spin-flip transi- 
tions between the unidirectional edge states can create 
nuclear spin polarization locally at the boundary of the 
2D sample^ - — Recently, we have predicted a similar 
phenomenon for a 2D TI, namely that embedded fixed 
spins such as the nuclear spins in a 2D TI can polar- 
ize locally at the boundary due to a current through the 
HESsj££ Interestingly, the 2D TI with localized spins re- 
mains ballistic , 57 ' 58 except if additional spin-flip mech- 
anisms for the localized spins are present^ 7 - However, 
combining localized spins and Rashba spin-orbit coupling 
in the 2D TI can produce a conductance change^ - — In 
the previous works; 5 - 7 - - — the interaction between the fixed 
spins embedded into the 2D TI and the HESs were mod- 
elled phcnomenologically. In contrast, here we pay spe- 
cial attention to the detailed forms of the HF interactions 
within a 2D TI. 

In this paper, we find the different HF interactions 
within the BHZ model for a HgTe QW. To this end, 
we take into account both the S- and P-like states of 
the BHZ model, which couple differently to the nuclear 
spins. We show that all the HF Hamiltonians couple 
the time-reversed blocks of the BHZ model. However, 
only HF interactions relevant for P-like states couple 
states within a time-reversed block as illustrated in FigfT] 
Moreover, we estimate the different HF coupling con- 
stants. The derived Hamiltonians arc general in the sense 
that they can be used to find the HF interactions for 
any kind of nanostructure in a HgTe QW, e.g. quantum 
dots^ ring structures^ quantum point contacts^ or 
hole structures^ 5 - As an illustrative example, we find the 
HF interactions for a pair of HESs. Remarkably, the intra 
HES transitions coupled to all the nuclear spin compo- 
nents perpendicular to the propagation direction of the 
HESs. This kind of coupling is unusual compared to e.g. 
an ordinary Heiscnbcrg model. Interestingly, the details 
of the HF interactions depend on the spacial direction of 
the boundary at which the HESs propagate. 

The paper is structured as follows: First the HF inter- 
actions and the BHZ model are outlined in sees. Inland 
IIIII Then the HF interactions are found within the BHZ 
model for the simplest case of a 2D QW fsec lIV[) . From 
this, we derive the HF interactions for a given nanos- 
tructure in sec. [V] Finally, the HF interactions for the 
HESs arc found and discussed (sec. IVI[) . Appendices lAllEl 
provide various details for completeness. 



II. THE HYPERFINE INTERACTIONS 

The HF interaction between an electron at position r 
with spin S = (S x ,S y ,S z ) and the nuclear spin I n = 
{Ix,n, Iy,n, Iz,n) of the lattice atom at R n can be derived 
from the Dirac equatio n 31 ' 47 to be (in SI units) 



, n MO 87r c/ s c T 

1 = 4^ Y leljn ^ ' ' 
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where /i™ is the Fermi contact interaction^, hJ]: is the 
dipole-dipole like coupling between the electrons spin and 
the nuclear spin, and h^ is the coupling of the electrons 
orbital momentum L n = r„ x p and the nuclear spin. 
Here r„ = r — R„ is the electrons position relative to the 
nth nucleus, r„ = |r„|, e„ = r„/r„ and u is the vacuum 
permeability. The gyromagnctic ratios of the electron 
7 e and the nth nuclear spin yj n of the isotope j are, 
respectively, given by (e > 0) 7 e = <? e us/fi,, where [Lb = 
eh/(2m e ) is the Bohr magneton and g e ~ 2 the electron 
g-factor; and jj n = gj n fi>N /h~, where fiN = eh/(2m p ) = 
/is/1836 is the nuclear magneton and gj n the g-factor of 
the jth isotope. Here m e and m p are the bare electron 
and proton mass, respectively^ 7 - Moreover, r c is a length 
scale related to the finite size of the nucleus and therefore 
much smaller than all other length scales in the system. 
It can be found to be^L r c = Ze 2 /(2mc 2 ) ~ Z x 1.5fm, 
where Z is the number of protons in the nucleus^ 8 - Thus, 
the total HF interaction between an electron and all the 
nuclear spins in the lattice is 



H 



HF 



Phf.i + Phf,2 + Phf,3 



(2) 



where only those lattice points R n with a non-zero nu- 
clear spin is included in the sum. 

Not every atom in a HgTe crystal has a non-zero nu- 
clear spin in contrast to e.g. GaAs. The amount of stable 
isotopes with a non-zero spin in Hg and Te are about^ 

17% of 199 Hg (spin-1/2), 13% of 201 Hg (spin-3/2), 
1% of 123 Tc (spin-1/2), 7% of 125 Tc (spin-1/2). (3) 

Hence, about 19% of all the atoms in HgTe have a non- 
zero nuclear spin. By isotope selection processes, this 
number can be varied somewhat experimentally. 

The contact interaction Phf,i is the only important 
HF interaction for S'-like states due to their spherical 
symmetry around the atomic core. On the other hand, 
P-like states vanish at the atomic core and therefore 
the contact interaction does not affect electrons in those 
states. In contrast, the two other terms Phf,2 arid Phf.3 



can indeed play a role for P-like states such as heavy- 
holes^ Moreover, Fischer et al£l found the atomic HF 
coupling constants to be about one order of magnitude 
lower for P-like compared to S-like states in GaAs. 



III. THE BERNEVIG-HUGHES-ZHANG (BHZ) 
MODEL 



Bernevig, Hughes and Zhang 2 ^ constructed a simple 
model describing the basic physics of a HgTe QW. The 
effective 4x4 BHZ Hamiltonian is derived using k • p 
methods^— and valid for k = (k x ,k y ) close to the T 
point, i.e. close to k = (0,0). The basis states of the 
model are the two Kramer pairs \E±) and \H±). Details 
on the derivation of the BHZ model arc found in Rcfs. 1, 
[H, and[ll For a 2D QW the BHZ Hamiltonian is 



fto = y>k#o(k)c k) 



k 



where c 



(4 



,,t 



k,_E+' l 'k,ff4 



-k,_E-' c k,ff-, 



(4a) 



is a vector of 



creation operators and 



ffo(k) 



( h(k) 
^ h*(-k 

with being a zero 2x2 matrix and 

h(k) = ( £k + ^ Ik A(ykx + iky ^ 
1 ' \ A(k x - ik y ) e k - M k 



(4b) 



(4c) 



Here e k = -Dk 2 , M k = M - Bk 2 , and k = Jk 2 + k 2 

have been introduced^ The parameters A, B, D and 
Mq depend on the QW geometry^^ Importantly, vary- 
ing the QW width changes the sign of Mq, which in turn 
makes the system go from a non-topological to a topo- 
logical state with HESs^ 

The hamiltonian (|4ap a priori has periodic boundary 
conditions and thereby does not contain any edges. By 
introducing boundaries into the model, it is possible to 
derive explicitly the HESs in the TI state of the QWH& 
This will be discussed further in sec. IVIAI 

Within the envelope function approximation^— the 
states of the BHZ model are 

\E+) = f E r 6 (z)\T 6 ,+l/2) +f Er8 (z)\r 8 ,+l/2), (5a) 

|#+> = M0|r 8 ,+3/2>, (5b) 

\E-) = f ET6 (z)\T 6 ,-l/2) + f Ers (z)\T s ,-l/2), (5c) 
\H-) = f H (z)\r 8 ,-3/2), (5d) 

where fi(z) are the transverse envelope functions in the 
z-direction perpendicular to the 2D QW and \Ti,m,j) 
are the lattice periodic functions^ at k = for the Ti 
band with projection rrij of the total angular momentum, 
J = L+S, on the z-axis. Here S is the electron spin and L 
is the orbital angular momentum (sec Appendix [Al . The 



time-reversal operator connects states within a Kramer 
pair (Q\E±) = t\Et) and @\H±) = t|#T)), and the 
two blocks in Po(k) (|4b|) are related by time-reversal. 
Here we choose phase conventions of the envelope func- 
tions such that time-reversed partners have equal enve- 
lope functions. Moreover, f E r 6 and fu are chosen real, 
whereas f E r 8 is chosen purely imaginary. (Appendix [A] 
gives more details on the envelope functions and the lat- 
tice periodic functions.) 

The states |£±) are seen to be mixtures of the .S-hkc Tq 
band and the P-like T$ band with rrij = ±1/2, whereas 
\H±) consist only of the P-like Tg band with rrij = ±3/2. 
Hence, the states have a definite total angular momentum 
projection, 



J z \E±)=±-h\E±) 



and J Z \H±) =±-h\H±), (6) 



but |P±) are not eigenstates of J 2 . 

The HF interactions can only induce transitions be- 
tween states with a difference of angular momentum pro- 
jection of one unit: rrij — rriji = ±1. Therefore, we can 
already at this point see that only particular combina- 
tions of the BHZ states can be connected by HF inter- 
actions as seen in Fig. [TJ Furthermore, it is evident that 
HF interactions relevant for both S- and P-like states 
need to be included to have a full description of the HF 
interactions in a HgTe TI. 

The real-space basis functions of Ho (|4"a|) for the 2D 
QW with periodic boundary conditions are 



Vk,E±( r ) = 



yJL x L y 



f „i{kxX+kyy) 



(7a) 



fErA z ) u v. +i( r ) + /sr 8 ( z ) u r a .±i( r ) 



¥>k,ff±( r ) 



ET 
V 



y L> x Ly 



- e i^ x+ k y y) fH{z) u {r ), (7b) 



where r = (x, y, z), L x (L y ) is the QW length in the 
x (y) direction and ur i ,m j {^) = (r|Ti, 7tj,j) are the real- 
space lattice periodic functions at k = 0. Moreover, we 
have included the atomic volume^ v a explicitly here as 
it is often done for HF related calculations 4i~— It de- 
pends on the choice of the individual normalization of 
the envelope functions and the lattice periodic functions, 
respectively, if v a should be included explicitly^ as dis- 
cussed in Appendix |B| 



IV. HYPERFINE INTERACTIONS WITHIN 
THE BHZ MODEL 



Next, we find the HF interactions within the BHZ 
model by using the states (O for a 2D QW with peri- 
odic boundary conditions. As we shall see, these results 
are useful, since they allow us to find the HF interactions 
for any nanostructure created in a HgTe QW (sec. |V|) . 



Outline of the way to find the hyperfine 
interaction matrix elements 



The HF interactions ([1} are local in space on the 
atomic scale, so the important part of the wavefunc- 
tion with respect to the HF interactions is the behav- 
ior around the nucleus. Hence, in the envelope function 
approximation, it is the rapidly varying lattice periodic 
functions ttr^m ■ (r) that play the central role, whereas the 
slowly varying envelope functions only are multiplicative 
factors at the atomic nucleus, as we shall sec below. 

We set out to find the HF interactions 

%IF,i = ^ E (¥^«r|-ffHF,i|Vk' V 'r')4«r C k't;'T' ( 8 ) 



k,k' 



=± 



for i — 1,2,3 in the basis ([7]), i.e. for </9k,i>±( r ) 
with v = E,H. We begin by describing the general 
way that we find the HF interaction matrix elements 
( i PkvT\H'HF,i\ i Pk'v'T')- To this end, the integration over 
the entire system volume V is rewritten as a sum of in- 
tegrals over each unit cell m of volume Vu' c , i.e. 



dv{- 



E 



d P {- 



(9) 



This should be understood in the following way: Every 
space point r can be reached by first a Bravais lattice 
vector Tt m = (X rn ,y rn ,Z m ) and then a vector p within 
the m th unit cell, i.e. r = TZ m + p. The superscript (m) 

on the unit cell volume w,™ indicates that the integral 
is over the mth unit cell. Thus, the matrix element is 

('PkvT\HnF,i\ l Pk'v'T') (10) 

= EE [, dpvi VT (K m +p)h^ w , v , T ,{ii m + p). 

Here one sum is over all unit cells 7L m , whereas the other 
sum is only over those atoms at position R„ with a non- 
zero nuclear spiniSi To proceed, we take v = v' = H as 
an illustrative example and obtain 

LxLy n TL m 



X \fH( Z m)\ / Jp Ur 8 ,r3/2(P) ft >r 8 ,r'3/2(p), 

where we have used the slow variation of the envelope 
functions on the atomic scale, /h(Pz + Z m ) ~ /jj(.£ m ), 
and the lattice periodicity of the lattice periodic func- 
tions, e.g. u rs _ r 3(72. m + p) = u rgT |(p) for all H m . Here 



n„ 



(X m , y m ), and the integral over p is over the m 



th 



unit cell, whereas hf is for the n th nuclei. For a specific 
nuclear spin n, we now include only the integral over that 
particular unit cell containing the n th nuclear spin, since 
the HF interactions are local in space. In other words, if 



the nuclei spin n is not inside the integration volume of 
the unit cell m, then the contribution is neglected^ i.e. 

(<PkH T \HBF,i\<PvHT>) = -^Ve'( k '- k )^ (12) 

y n 

x \f H (Z n )\ 2 / dpu* r8T3/2 {p)h?u raT , 3/2 (p), 



T e ,m.j 


(13a) 


a Te *?:,m, (r+ g ) - «H g *r 6 S mj ('- 


■!)] 


i 


r 8 ,m 3 - 


(13b) 


«Tc*?:^( r+ t) +aH ^"%( r - 


■!)] 


i 



where the unit cell integral now is independent of the 
unit cell position TL n . The sum is only over the lattice 
nuclei at R„ with a non-zero nuclear spin. Therefore, the 
system does not have discrete translational symmetry, so 
the sum cannot simply be made into an integral. Hence, 
the matrix elements are not diagonal in k due to the 
nuclear spins at random lattice points. 

In order to proceed, we need to evaluate the integral of 
the lattice periodic function over the unit cell in Eq. (fT2")l . 
To this end, the symmetry of the lattice periodic func- 
tions arc important: The contact interaction 7Jhf,i is 
zero for P-like states, since they vanish on the atomic 
center, while matrix elements of i?HF,« for % = 2, 3 vanish 
for 5-like states due to their spherical symmetry. Here, 
we approximate the lattice periodic functions by a Linear 
Combination of Atomic Orbitals (LCAO) as^o 

u r 6 .m» = N 1 



u r s . m ,(r) =N } 



where ^y c m . and ^ r g m . are atomic-like wave functions 
centered on the Te and Hg atom, respectively, and r 
is only within a single two-atomic primitive unit cell of 
HgTe centered at r = 0. The atomic wave functions in- 
herit the symmetry of the ban d 44 ' 80 as indicated by the 
index r,-, m,j. The atoms are connected by the vector d, 
and the constants Nr i , m - are determined by the lattice 
periodic function normalization J dr\uT i ,mA r )\ 2 = 2, 
see e.g. Eq. (|B4j) . The electron sharing within the unit cell 
is described by a Te , H •., which fulfill IcktcP + I^HkI 2 = ^^ 

The LCAO approach (fl"3"|) now facilitates evalua- 
tion of the unit cell integral in the matrix elements 
(ipkvT\H'HF,i\ l P'k'v'T') ■ Consider e.g. the unit cell inte- 
gral in Eq. (fT2"j) for a non-zero spin on the n th Hg nucleus 
located on p = d/2, i.e. 

/ d P U *r 8 .r3/2(P) h ? U T 8 ,T'3/2{p) (14) 

- N r 8 ,T3/2 N r8,T'3/2\u Hg \ 2 

x I dp[^l T3/2 (p-d/2)]*h^ T , 3/2 (p-d/2), 

where only the important contribution of the atomic wave 
functions centered on the Hg atom is included. In other 



words, integrals involving atomic wave functions centered 
on different atoms are neglected. Fischer et al^ esti- 
mated that these non-local contributions are two to three 
orders of magnitude smaller for GaAs - even for the long- 
ranged potentials in h^ 3 in Eqs. ([rbllc[) . 

Thus, we have now outlined how to find the ma- 
trix elements ((pkHrlH-nF^lfk 1 Ht') for all three kinds of 
HF interactions (flj. The matrix elements of the types 
(VkEr|-H r HF,j|</'k'ET') and ((P^etIH-hfA^'Ht') follow the 
same lines as above. The essential ingredients arc the lo- 
cality of the HF interactions, the periodicity of «r„m (r) 
and the slowly varying envelope functions. Next, we find 
the three HF interactions (JXJ) within the BHZ model. 



B. The contact HF interaction for S'-like states 



Now we find the contact HF interaction Phf.i Eq. fTaj) 
within the BHZ basis (J7]). We begin by noting that 

^Ht\Hb_F,iW^'Ht') = and (^kErl^HF.ll^k'Hr') = 0, 

since the contact interaction is only non-zero at the 
atomic center (ft™ oc S(r — R„)), where the r§ P-like 
atomic orbitals vanish. Hence, only the Tq S'-like part 
of the fk'Er'i^) states leads to non-zero matrix elements 
of -Hhf,i- Using the approach in Sec. IIVAI to find the 
matrix elements, we get 

(<PkEr\H SF ,l\'PkfEr') = -^ Ve'( k '- k >' K - (15) 

l^x^y — 
a n 

X \f E r 6 ( Z n)\ 2 / rfP'"r 6 ,rl/2(P) /l l U r6,r'l/2(P) 
Jv ua 

for t,t' = ±. The T 6 states ur 6 .±i/2( r ) simply factor- 
ize into a spin and an orbital part as ur 6 .+(-)i/2( r ) = 
ur 6 (r)| t (I)), SCG e.g. Eq. (|A2|l . Using this and the ex- 
plicit form of the contact interaction ft™ in Eg. ([Taj) , we 
readily obtain 



(<PkET | PhF, 1 1 <Pk' Et' 



X S,j n (Zn) 



±j x Lj 



Ve. 



,i(k'-k)--R.„ 



e 

1 . „ 1 



(16) 



' cj ^-Z,n^T,T' ~T n t' tTh^T^—T 1 





199 Hg 


201 Hg 


123 Tc 


125 Te 


A f°r c i^ eV ] 


4.1 


-1.5 


-49 


-59 


A R°r c t/ ieV ] 


0.6 


-0.2 


-6.0 


-7.2 



TABLE I. Estimates of the atomic contact HF couplings 
^ Atomic Eq.pgji and the atomic P-like HF couplings ^ A '°" lic 
Eq. (|26[) in HgTe for the naturally present isotopes with non- 
zero spin, see ((3|. The HF couplings for S'-like states are seen 
to be about one order of magnitude larger than for P-like 
states. The sign of the HF couplings stem from the sign of 
the nuclear g-factors. See Appendix [D] for details of these 
estimates. 



tion, since it can be given in terms of the atomic or- 



bital *f" 



by using Eq.flU as A At ? mic ex |u rr (R„)| 2 



M 



r 6 ,i/2i 



i*f; 



(0)| 5 



i.e. 



A Atomic 
A S,j n 



only depends on 



the nuclear isotope type j n at site n. Moreover, at the 
present level of approximation, we can freely replace the 
Bravais lattice vector TL n by the actual position of a nu- 
clear spin within the n th unit cell in the envelope func- 
tions in Eq. (|16[) due to their slow variation. Finally, we 
arrive at the HF contact interaction in the BHZ basis as 



^HF, 



e i(k'-k)--K nX 

12 12 —T-T C k#HF,lC k , 

n k,k' 



L X Ly 



(19) 



where c 



J 



J 



J 



J 



,C k,_B+' C k,ff+' C k,E-' C k,ff-J and— 



Hhf,i — WZ A S,j Tl 



1 

2ft" 



(Zn) 



h,n J_,„ 



/+,„ -J,, n 





(20) 



The sum is only over non-zero nuclear spins. Therefore, 
it is now clear that the contact HF interaction contains 
elements ex I± t n, which connect the time- reversed blocks 
in the BHZ hamiltonian (|4bl) . Moreover, as illustrated 
in Figfl] only the |P±) states are connected by Phf.i, 
since only these states contain a S'-like symmetry part. In 
in table HI estimates of the atomic contact HF couplings 
^Atomic are gj ven f or t ne stable isotopes of HgTe with 

non-zero nuclear spin (see Appendix ID1 for details). 



where i±, n = I x ,n^Hy,n are the raising and lowering nu- 
clear spin operators. In analogue to the case of a quan- 
tum dot^^ we here introduce the position dependent 
contact HF coupling as^ 2 - 



A StJn (Z n ) = v a \f ET6 (Z n )fA%?™, 
which includes the atomic contact HF coupling 



A 



Atomic 

S,jn = 



2/lQ 



9e^Bg 3n fJ-N\u rf ,(Rn) 



(17) 



(18) 



for the nuclear spin at site n of isotope j. Here A s (Z n ) 
depends on the real space position of the nuclear spin. In 
contrast, Ag*"™ 110 does not depend on the nuclear posi- 



C. The HF interactions for P-like states 

Next, we find the HF interactions within the BHZ basis 
d?]) for i?HF,2 and #hf,3 Eqs. (|lbllc[) . which are relevant 
for the P-like states. 

To begin with, we argue that the Tq S-likc states - 
part of the E± states - do not contribute to the matrix 
elements ((pkEr\HHF,i\'Pk'ET') and (93kEr|-HHF,i|</>k'.Hv) 
for i = 2,3. (In contrast, the Tg P-like part of E± do 
contribute to these elements as will be shown below.) To 
understand this, the HF matrix elements are written in 
terms of the unit cell integrals over the atomic-like wave 
functions as outlined in Sec. lIVAl Firstly, for the dipolc- 



dipolc like HF interaction (jlb|) . we have 
/ dp[^°( P Td/2)]*h^ffipTd/2) = (21) 

due to the rotational symmetry of the S'-like orbitals 
around the atomic core.— Secondly, we have 

/ dp[^ T ;(p T d/2)]*^<2( / ,Td/2) = (22) 

Jv uc J 

due to opposite parities of the S- and P-like orbitals.— 
The same matrix elements containing h^ instead of h^ 
are also zero, because the .S'-like states have zero orbital 
momentum, i.e. L n * r ^ / m c (p) = 0. 

Therefore, only P-like states contribute, so we are now 
left with (see Sec. HVA"|) 
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L X Ly 
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(23a) 



x\f E r 8 (Z n )\ 2 / dpUr 8)T i /2 (p)^>r 8 ,r'i/2(p) ! 
{<PkHT\Hm?, i \<p k ,HT>) = j^r E e t(k '- k) ' TC - (23b) 

y n 

x\f H (Z n )\ 2 / dp Ur 8 ,T3/ 2 (p)fc"«r 8 ,T'3/2(p), 

(^kSrl^HF.ikk'ffr') = y^r- V e*( k '- k )^ (23c) 

x /Br 8 ( z «)/ff( z «) / d P«r 8 ,rl/2(P) /l ™ U r s ,T'3/2(p) 



and (^kHrl^HF.il^k'Br') = (^k'Br'I^HF^I^kffr)*, 

where i = 2,3 and r, T 1 = ±. Using the LCAO 
approach Eq. (fT3"|) . the unit cell integrals over the lat- 
tice periodic functions now become integrals over the 
atomic- like wave functions as in Eq. (TTJJ. We write 
the atomic wave functions as a product of a radial 
part i? Hg / Tc (r) and an angular part Yr 8 , m .(#, 4>), i.e. 

*r 8 S/ J>) = R Ks/Tc (r)Yr 8 . m] (d,cf,), using spherical co- 
ordinates (r, 6, <f>) with the nucleus in the center. Since 
the integrals are over the two-atomic unit cell volume, 
they do not a priori factorize into a product of radial 
and angular integrals. However, due to the 1/r 3 depen- 
dence of hf (i = 2,3), the important part of the unit cell 
integrals are numerically within one or two Bohr radii ao 
from the atomic core, which is certainly within the unit 
cell volume. Therefore, it is a good approximation to 
write the unit cell integrals (e.g. Eq. (|T4")0 as 



Jv uc J 



(24) 



drr 2 j d^ j d9sm(8){^l m (r)yh^l ml (r), 



,Hg f_\i*z.ra,T,Hf> 

>o Jo Jo 

where the specific choice of r max ^ &o is n ot important 



for the numerical value of the integral^ Therefore, we 
are now left with an essentially atomic physics problem, 
where the integral separates into a product of a radial 
and an angular part. The radial part is 



1 



drr 2 \R K (r)\' 



r*(l + *y 



(25) 



which is the same for all the matrix elements of h^ and 
/13 and only depends on the type of atom k =Hg, Te. 
Due to the smallness of the nuclear length scale r c , it is 
not significant for the magnitude of (l/r 3 )£i££ Using the 
radial integral (|25p . we introduce the atomic P-like HF 
coupling for isotope j (at site n) as 
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Atom 
P,3n 



Mo 
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ff e ^sff ? „MAr(^r 8 ) 2 |aj 



1 



(26) 



which are estimated to be about one order of magnitude 
smaller than the atomic contact HF couplings ^4^*,™"° 
(TT51) . sec table HI Here it makes sense to have a com- 
mon atomic HF coupling for the dipole-dipole like cou- 
pling h% and the orbital to nuclear-spin coupling /ig, 
since the normalization constants for the LCAO lat- 
tice functions (|13[) are numerically approximately equal, 
iVr 8 .3/2 — Ar 8jl / 2 = Ar s , as discussed in Appendix |Pl££ 
Moreover, we also use that Nr iyln - are independent of the 
sign of rrij , see Eq. (|D6|) . Calculating the angular integrals 
as discussed in Appendix [C] the matrix elements (|23l) for 
the dipole-dipole like HF interaction Phf,2 become 
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(27c) 



where we introduce the position dependent P-like HF 
couplings as 



aHH _ \f cpr \ 124 Atomic 
aEE if ly \|2 /(Atomic 
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(28a) 
(28b) 
(28c) 



and APj = [Apt ]*. In comparison, for the con- 
tact HF interaction only a single position dependent 
HF coupling was introduced in Eq. (fl"7|) . Here the ex- 
plicit dependence on the position Z n of the nuclear spin 
has been suppressed in the notation for simplicity, i.e. 



A v p- (Z n ) = A v p- . Similarly, the matrix elements for 
the HF interaction Phf.3 between the electronic orbital 
momentum and the nuclear spins become 

(Vk^rl^HF.slVk'Sr') = Y\ -i-e^'-^ 7 ^ (29a) 

L x L y 
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(29b) 
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(29c) 



It is noteworthy that the heavy hole like states H± only 
couple diagonally (r = t') or Ising-like in (|27bl29b[) in 



agreement with Ref.[4J|. Physically, this is because the 
H± states have a difference of total angular momentum 
projection larger than one, \rrij — rriji\ > 1. Moreover, 
the coupling between the states E± in Eqs. (|27al29a[) is 
essentially like the coupling between the light hole states 
|r§,±l/2), since the S-like states do not contribute to 
the matrix elements of Phf.2 and Phf,3- For these ma- 
trix elements between the E± states, the off-diagonal 
elements (r = — r') are a factor of 2 larger than the di- 
agonal elements {t = r') in accordance with Ref . j48j . 

Using the matrix elements in Eqs. (|27l29[) . we now fi- 
nally arrive at the HF interactions relevant for the P-like 
states in the basis ([7} as 
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for i = 2,3, where 
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(31a) 



and 



PhF,3 — — 5i?HF.2 



(31b) 



V. HYPERFINE INTERACTIONS FOR A 

NANOSTRUCTURE IN A HgTe QUANTUM 

WELL 



such that the total HF interaction for the P-like states 
becomes 
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HF,P 
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HF.2 
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HF,3 



AH 



HF.2- 



(32) 



Just as the contact HF interaction (|20p . the P-like HF 
interaction connects the time-reversed blocks by connect- 
ing the E± states. Moreover, the P-like HF interac- 
tion connects the states within the time-reversed blocks 
(e.g. E+ and H+) in contrast to the contact HF inter- 
action, see Fig. [1] 



Interestingly, the sign of the dipolc-dipolc like HF in- 
teraction (|31ap is opposite to the contact HF interaction 
(|20|) and to the orbital to nuclear-spin coupling (I31bl) . 
However, since the elements of Phf,3 are larger than 
those of Phf,2 in absolute value, the total HF interaction 
for the P-like states (|3"2"j) end up having the same sign as 
the contact HF interaction. 



Now, we show how the HF interactions for any nanos- 
tructure in a HgTe QW can be derived from our results in 
Eqs. (|2QI31[) for a HgTe QW with periodic boundary con- 
ditions. Examples of such structures are quantum dots^ 
mesoscopic rings^ point contacts^ and anti-dots^ 



For a given nanostructurc, the envelope wave func- 
tions are needed to find its HF interactions within 
the BHZ framework. Utilizing the Peierls substitution 
(k x ,k y ) = —i{d x ,dy) in the BHZ hamiltonian (|4b|). 
the envelope functions <I> (rj_) can be found by solv- 
ing [H {-id x , -id y ) + V(r±)]$ ri (r 1 _) = E ri ^ v {r ± ), where 
V(r±) is the potential confining the nanostructure,— 
r_L = (x, y) and r\ is a collection of quantum numbers 
to be specified for a concrete nanostructurei 62 ' 63 ' 65 ' 74 i 88 
The envelope function is given by 8 ^ 



%(r±) = 



( <t> v ,E+( V ±) \ 
4> v ,E-( r ±-) 



(33) 



such that the entire wave function including the lattice 



periodic functions is 



iMr) = V^ E ^,c( r ^)( r K>: (34) 

C=B±,ff± 



a semi-infinite half plane restricted to x > and periodic 
boundary condition in the y-direction, k y is still a good 
quantum number. The HESs envelope functions running 
along the y-direction become^ (see FigJ2|) 



where (r\E±) = f E r 6 (z)u r ^ ± i{r) + f E r 8 (z)u rs>± i(r) 
and (r\H±) = f H (z)u r ± 3 (r). For instance, the 2D QW 
with periodic boundary conditions simply has the en- 
velope functions e l ( k * x+k y y) (I, 0, 0, Q) T j ^TJ7 y etc., see 
Eq.©. In situations with time-reversal symmetry, the 
states appear in Kramers pairs of equal energy. The BHZ 
model is constructed such that Kramers pairs appear as 
equal energy solutions of the upper and lower 2x2 blocks 
in H a (SB. Thus, 



'o ii^u|i - j. mis, cp„ E _ = </>„ H _ = for one of the two 
states in a Kramers pair and vice versa, which simpli- 
fies the algebraic burden of finding $ , 62 ' 63 i 65 ' 74 The two 
states in a Kramers pair are sometimes referred to as spin 
up and down, since the upper (lower) block only consists 
of orbital states with positive (negative) total angular 
momentum projection, see Eq.©^ Here the index for 
the Kramers pair is included in r\. 

To find the HF interactions (JlJ for a given nanos- 
tructure with envelope wave function $ (rj_) (|3U)) . the 
same idea of separation of length scales as in Sec. IIVI is 
used: The HF interactions act on the atomic length scale 
such that the slowly varying envelope functions only be- 
come multiplicative factors in the HF interactions for the 
nanostructurc. Thus, we find 






(x,y) 



l 



e lkyV i)l (x)xl for s = u,d (36) 



with the energy dispersions E% = Eo + hvQk y and Ef = 
Eq — Hvaky, respectively. The dispersions are exactly 
linear in the semi-infinite half plane model used here^ 
Both the velocity v = -\AB 2 - D 2 \A\/(hB) and E = 
—MqD/B are positive for realistic parameters*^ and L y 
is the length of the edge. The spinor parts of the HESs 
are independent of k y and given by 
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(37) 



with the normalization factor n = y/(B 
transverse part of the HESs are tffc (x) 

t)i(x) = g_ k (x), where 
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(38a) 
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n rj.r}' 



< ^ > r;'('^-n-L) c 77 c 77': (35) with the k y - dependence inside Ai and A 2 as 



where i?HF,i are the 4x4 matrices found in Eqs. (|2QI31|32j) 
for the contact (i = 1) and P-like (i = 2, 3, P) HF inter- 
actions, respectively. The sum is only over the atomic 
sites n with a non-zero nuclear spin. The Bravais lattice 
vector TZ n± = (X n ,y n ) pointing to the unit cell contain- 
ing the n th nuclear spin can freely be interchanged by the 
atomic position R„ of the nuclear spin due to the slow 
variation of the envelope functions on the atomic scaled 
Thus, we have now given the general form of the HF in- 
teractions in the BHZ model for a given nanostructurc. 
Below we give an example. 
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Here A 2 1 determines the decay length scale of the HES 
into the bulk and 
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VI. HYPERFINE INTERACTIONS FOR A PAIR 
OF HELICAL EDGE STATES 

Next, we deal with the HF interactions for a pair of 
HESs in a HgTc 2D TI QW. 



A. The helical edge states along the y-axis 

To find the HF interactions, we first give the envelope 
wave functions for a pair of HESs. These can be found by 
introducing a boundary in the BHZ model and requiring 
that the envelope functions vanish at the boundary^ For 



The HESs only exist in the topological regime of the BHZ 
model where M$/B > 0. The explicit forms above were 
derived under the assumption < M /B < A 2 /(4_B 2 ), 
where Ai 2 are purely reaL 74 ' 92 This is the relevant regime 
for the realistic parameters 1 ^ for 2D TI in a HgTc QW 
of width 6lA or 70A. 

Using the time-reversal properties of the basis states of 
the BHZ model (as discussed in Appendix [A"f . it is seen 
explicitly that \^ fc (x, y) and ^ y .-k ( x >y) constitute a 
Kramers pair, since they are connected by the time- 
reversal operator O as Q^ y k (x, y) = —^y _ k (x, y) and 



&Kk y (^y) 



\T/" 

y, 



(x,y). Oftenia*?, « fe J is 



referred to as the spin-up (spin-down) edge state, since 



it only consists of states with positive (negative) total 
angular momentum projection, see Eq. p7[) . 



where 



-k v d 



are the creation (annihilation) opcra- 
36|) . We have emphasized 



tors for the HESs ** fc (x, y) 

(v) 
in the notation that n^F i is for HESs along the y-axis. 

The product of the transverse parts of the HESs at the 

nuclear spin n is introduced as 



B. Hyperfine interactions for the helical edge 
states along the y-axis 



The HF interactions are now readily found by inserting 
the envelope HESs along the y-axis (|36p into the general 
HF interaction formula (|35|) for any structure in a HgTe 
QW. Using Eq. ([2T))) the contact HF interaction becomes 

w& , =i^ E E —t — *>*&) ( 39 ) 



L HF.l 



2h 2B 



\( X n) T t _ a(*») 7 ^ r 

iV fc y ,fe^- , 2,™ C fe y « C fc;M Iy -k v -k' y 1 z,nC kyd C k ,^ d 

\ {Xn) T r f r -At*") r, r 1 " r 

Ii -k v ,-k'„ 1 -,™ c k v u c k'd li --k v ,k' 1 + ,n c k v d c k'u 



A Si=ift)ft;ft). 



(40) 



/(w) 



and includes the only dependence of X n in "H^p x . Here 
we see that the contact HF interactions can produce 
transitions between the HESs \fr" k (x, y) and ^ k , (x, y) 
at the expense of a change in a nuclear spin state. In 
particular, elastic transitions within the Kramers pair 
^y k ( x > v) an d ^t -k ( x j y) arc possible. This is just as 
if the HESs were spin- 1/2 as used e.g. in Refs.l57l.l58l. and 
I6G . Hence, from the point of view of the electrons in the 
HESs the time reversal symmetry is broken. Of course, 
the composed system of electrons and nuclear spins is 
time-reversal invariant, since any system can be made 
time-reversal invariant by expanding it sufficiently. 93 

(v) 
The HF interactions due to the P-like states, "Hjjp p = 



Wg$ 2 + ^hf 3' i s similarly found by inserting ty q k l[56]) 
into Eg . (1551 and using Phf.p l[52"]). i.e. 
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where the dependence on Z n is inside the HF cou- 
plings ApJ Eg. ([28)) . Here we used the rewrit- 
ings ^I+.Jifi - I-.nAp'Zj = -24,Jm(A|fJ and 



i[I+,„^ B 



PjV, 
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on the fact that Af,*j 



-.nAfVj = +2I x , n Im(Af H 



ilm(Af H 



fljj, which build 
Pj n j due to the phase 
conventions of fEr 8 as purely imaginary and fjj as real. 
This HF interaction also permits transitions between the 
two HESs - especially within the Kramers pair - just as 
the contact HF interaction (|M|) . The terms oc cLc fe/? in 
the HF interaction (|4ip affect transitions within a single 
HES. These are more unusual than their counterparts in 
the contact HF interaction (|5§)) . since they do not only 
contain terms involving I z , n , but also I x ,n- Hence, the 
HF interaction (|4"Tj) due to the P-like states has terms 
like a I X S Z coupling, which are not present in e.g. a 
Heisenbcrg model. These terms oc I Xjn (y k c k , stem from 
the fact that the HF interactions due to the P-like states 
Eq. tptTj) couple the states |P±) and |P±) within a single 
time-reversed block of Hq. 

In order to shine more light on the form of the HF 



I 

interactions ([M| and (|4"Tj) , they are given in Appendix [E] 
in terms of non-diagonal edge state spin operators using 
the picture of spin-1/2 HESs. 



C. Hyperfine interactions for helical edge states 
along the s-axis: Curious differences 

The HF interactions presented above are for HESs run- 
ning along the y-axis. Now, we find various interesting 
differences in the HF interactions for HESs running along 
the x-axis. 

The HESs are found in the same way as in sec. IVI Al 
The only difference is that we consider the HESs localized 
near a boundary given by the x-axis instead of the y-axis, 
i.e. we study the semi-infinite half plane defined by y > 0. 
The HESs along the x-axis are given by 
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FIG. 2. The helical edge states along the j/-axis Eq. (|36p and 
along the a;-axis Eq. ()42[) are connected by k y -H- —k x . More- 
over, they differ by an imaginary unit i in the E spinor compo- 
nents [compare Eqs. (|37[) and (|43[) ]. This leads to interesting 
differences in their HF Hamiltonians (see the main text). For 
illustrative purposes, the HESs are drawn side by side even 
though equal energy HESs in fact are on top of each other. 
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i.e. the imaginary unit i does not appear in the E± com- 
ponents of the spinors as for the HESs along the y-axis, 
see Eq. (|37[) . This is the mathematical origin of the dif- 

I 



ferences between the HF interactions for the HESs in the 
two directions. These HESs also appear in Kramers pairs 



(*£*. and <_ fci 



and $>1 k 



(^x fc x ) i s referred to as 



spin-up (spin-down). The spin- up HES ^" k has nega- 
tive velocity such that E k = Eq — hvok Xl while the spin- 
down HES has positive velocity, i.e. E k = Eq + hvok x . 
Hence, the velocities of the HESs along the x and y axes 
have opposite signs, such that spin-^ always travels the 
same way along the boundary, see Fig. O Therefore, it is 
natural that k x has to be exchanged by — k y to connect 
the HESs in the two perpendicular directions. 

By inserting the HESs along the x-axis Eq. (|4"!Z]) into 
Eq. (|35|) . the contact HF interaction becomes 
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dtjy^gk'S^- Interestingly, the sign 
of the terms producing inter HES transitions is opposite 



39l) . This difference stems from the 



(v) 
to the one in "Hjjp i 

imaginary unit i in x y (J37J) , which is absent in x x lj43[) . 
Moreover, the sign of k x and k y is opposite in the A func- 
tions for ?4if,i © and H^-y d39]), respectively. This 
is natural in order to maintain the propagation direction 
of the HES-spin <r = u, d, see Fig. [2j 

The HF interaction due to the P-like states for the 
HESs at the cc-axis becomes 
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where the inter HES transition terms again have an 
opposite overall sign compared to H^p P (|4*T|) . An- 
other noteworthy difference is the exchange of the terms 

^," C fe y? C fc; ? in n< HF,P b y h-nA^k'^ in W HF,P> Lc - intra 
HES transitions are coupled to the nuclear spin operators 
perpendicular to the propagation direction. These differ- 
ences again stem from the imaginary unit (or the lack 
thereof) in the spinors. Furthermore, the signs of k x and 
k v are again interchanged in the A functions by compar- 
ing Wnl.p and "Hhf p- 



D. Position averaged hyperfine interactions 



In HgTc about 19% of the atoms have a non-zero nu- 
clear spin and these can be assumed to be randomly dis- 
tributed. In the HF interactions for the HESs the unit- 
cell position of every nuclear spin is included. This infor- 
mation is sample dependent and often valuable insights 
can be found without it. Therefore, we now consider the 
HF interactions averaged over the unit-cell position of the 
nuclear spins in analogue to impurity averaging.— To be 
specific, we focus here on the HF interactions (|39j) and 
(|4"TT) for the HESs along the y-axis. Mathematically, the 
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position averaged of some quantity JF is introduced as 



T . 



1 



A 



N, 



dlZo 



dno, N3 T(n 1 ,...,n Ns ), (46) 



where N s is the number of non-zero nuclear spins covered 
by the HESs. We only average over the positions TZ , n = 
(X n , Z n ) in the cross section area A of the HESs along the 
y-axis. Thereby we keep the positions y n , which break 
translational invariance along the edge and ultimately 
can lead to backscatteringi^2r— 

Now we study the position averaged HF Hamiltoni- 
ans. However, one can equally well position average at a 
later stage of a calculation, if it is physically relevant for 
a particular phenomenon, e.g. position averaging of the 
transition rates^ Using the normalization of fEr e in Ap- 
pendix [B] the position averaged contact HF interaction 
becomes 





199 Hg 


201 Hg 


123 Te 


125 Tc 


A* jn [MeV] 


1.1 

(77%) 


-0.38 

(75%) 


-12 

(73%) 


-14 

(73%) 


Af n [fieV] 


0.30 

(14%) 


-0.11 

(12%) 


-3.5 

(12%) 


-4.2 

(12%) 



TABLE II. Estimates of the effective HF couplings (J53J) for 
a pair of HESs. In parentheses we give the percentage of 
the HF coupling stemming from the HF hamiltonians due to 



P-like states, e.g. 4[(P 



D)/(15B)]A^ jn 



™/Aj 



Here we 



use the atomic HF couplings in table [J and the BHZ model 
parameters' 1 *^ B and D only for a 7uA thick QW. 



and k' y for D ^ 0, i.e. A ky>k > v ~ T] ky , ky , where 

D 2 \A 2 B+2(B 2 -D 2 )M ] , lx , , N 

^ = 1 [ 8 BM 2 (B 2 -D 2 ) l ^-K) 2 - (49) 



{y) _ i b-d sr^r 

Hn¥ > l ~4h 2B ^^ 



A Ate 



N 



(47) 



A-k y ,k'Iz : nC k C k , A._k v ,-k'Iz,nC k d c k>,d 



"■ky,-k'I-, n C k u C k , d A—k y ,k'I+,nC kd C k , u 



where the cross section area A is given in terms of the 
QW thickness W z and the HES width W x as A = W X W Z , 
such that J. dHo.i = 1. Here W x is on the order of a few 
decay lengths A^~ and the number of atoms covered by 
the HESs is N = (LA)/v a - We observe that the posi- 
tion dependent HF coupling A s ■ (Z n ) (fT7|) is replaced 

by a homogenous HF coupling Ag^^/N due to the po- 
sition averaging as in the case of quantum dots.— The 
position average of the product of transverse functions, 



-w x 



Afe,fe' = J " dX n g* k (X n )g k ,(X n ), is now independent of 
the positions X n . It can be well approximated by replac- 
ing W x by oo in the upper limit, which gives 



A ky,k' y 



, 2A 1 A 2 (A 1 +A 2 ) / 2A' 1 A 2 (A / 1 +A 2 ) 
" (Ai-A' 2 )2 



(Ai-A 2 ) 2 

1 
X 1 + Ai A 2 + Ai \ 1 + A 2 A 2 + A 2 



(48) 



1 



Hence, the lowest order expansion fulfills rj k jk i 



V-k u ,-k' and r;_ 



k y ,k 



r\ k — k i such that the position 



averaged contact HF interaction (|4T|) simplifies to 

in n 4 Atomic 

«S.i <4^ E E^-^%- ( 5 °) 



4ft 25 



Vk y ,k' y Iz,n(C kyU C k ^ u C kyd C ky d) 



Vk y -k'(I-,ncl u C k , d + I+, n c\ d C k , ) 



to lowest order in k y and k! . In this limit, n^p 1 therefore 
has uniaxial anisotropy. 

The position averaged HF interaction due to the P-like 
states in Eq. (j4T|) becomes 



u 



(y) 

HF,P 



IKS Z^ 2-^, 



15ft 




A Atomic 
c i{k' v -ky)y n - t± P,3n 

N 



(51) 



Iz,n[ C k v u C k'u C k v d C k'd\ 



»* n V "I" 5 



B ^ k v- k 'y 



*-~,n c k v u C k'd + I+,n C k v d C k'u 



where A^ and A^ depends on k y and k' respectively. It 
is evident that A.k v ,k' = A k ^ ky , since g ky in Eq. (|38a[) 
is real. Moreover, A ky , ky = 1 due to the normaliza- 
tion of g k . Furthermore, in the particle-hole sym- 
metric limit D = 0, we have g- k (x) = g k (x) such 



-(D=0) 



A_ 



-(D=0) 



Aa 



-(D=0) 



that Afc k i 

A- ky ,k' . Hence, the position averaged contact HF 

interaction (|47|) becomes isotropic in the particle-hole 
symmetric limit. Since the BHZ model is valid only close 
to the r point, we expand A k k i to lowest order in k y 



by using the expansion A ky . k ' ~ f]k v ,k' and the nor- 
malization conditions for fjj and /_Er 8 ( SCG Appendix 
|BJ) . Interestingly, the terms in Hj^ F P (HJ) coupling I XtU 
and c k c k , vanish in the position averaging, since 

Jh{z) is even and JeTs,{z) is od d 25 ' 72 such that Ap 1 ? oc 
J dZnfgY (Z n )fjj{Z n ) = 0. Furthermore, even in the 

particle-hole symmetric limit D = 0, H-j^p P is not 
isotropic in contrast to the contact HF interaction. 



The total position averaged HF interaction Ti^p 
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Wj![F j + "HfjF p in the small wave- vector limit is now 



found from Eqs. (|50|51[) to be 



H 



(y) 

HF 



1 V e l( - k 'y- k y )y " 



2h 



(52) 



.4 



r lky,k' y l-z,n\ c k y u C k' y u c k y d C k' y d) 



A± 



^- 7 lk y -k' y V-,nC kyU c k ^ d + I +,nC kyd C k ,^ u ) 



N 



where effective HF couplings were introduced as 



A] 



Ai 



,rB + SD #, (53a) 



.4 



Atomic 



45 Sj 
B -D 
B 



15 25 



^ 



1 4 

_ /(Atomic I /i Atomic 



V?,j„ 



15 



l P,Jn 



(53b) 



Hence, the total position averaged HF interaction "H^ F 
has uniaxial anisotropy. Estimates of A| and A^ are 
given in table [Til Remarkably, the part of the effective 
HF couplings due to the P-like states dominates for the 
coupling A?, but not for Aj- . One reason is that the 
HESs have their main contribution on the H states com- 
pared to the E states, since (Xy V - (-i0.36, -0.93, 0, 0) 
for a 70A thick QW^i Moreover, not only the H states 
are P-like states, but also partly the E states, see Eq.©. 

We remark that the position averaged HF interactions 
for the HESs along the x-axis Eqs. (|44l45|) follow along 
the same lines. The only difference in the total HF in- 
teraction in Eq. ([52|) is an opposite sign of the inter HES 
transition terms (apart from the replacements y n — > X n 
and k y — > k x ). 

For typical parameters^ the number of atoms covered 
by the HESs is about N ~ 10 7 per /xm edge, where about 
19% of these atoms have a non-zero nuclear spin. 



VII. DISCUSSION AND SUMMARY 

In summary, we have provided benchmark results 
within the BHZ model for the form of (i) the contact 
HF interaction in Eq.(j20|). (ii) the dipole-dipole like HF 
interaction in Eq. (|31a[) and (iii) the coupling of the elec- 
trons orbital momentum to the nuclear spin in Eq. ()31bl) . 
All the HF interactions couple the time-reversed blocks 
of the BHZ Hamiltonian - just as the Rashba spin-orbit 
coupling^ However, for HF interactions the coupling 
persists even to the T point (k = 0) and they break time- 
reversal symmetry from the electronic point of view. 

Furthermore, we estimated the atomic HF constants 
relevant for a HgTc QW, see table HI These estimates are 
generally smaller by an order of magnitude or so com- 
pared to similar estimates for GaAs by Fischer et a L 44 i 97 
This is natural, since heavier elements often have lower 
HF couplings due to their higher principal quantum num- 
ber of the outermost electron [see e.g. Eqs. (|D4ID5[) ]. 



From the HF Hamiltonians within the BHZ model, we 
derived a general formula (f3"5"j) for the HF interactions 
for any nanostructure in a HgTc QW. The input of this 
formula is the envelope function of the given structure. 
From this formula, we found the HF interactions for a 
pair of HESs. Interestingly, the HF Hamiltonians depend 
on the orientation of the boundary at which the HESs 
propagate: The sign of the terms creating inter HES tran- 
sitions is opposite for perpendicular boundaries. Further- 
more, the HF interactions due to the P-like states couple 
the intra HES transitions to both nuclear spin compo- 
nents perpendicular to the propagation direction of the 
HESs, see Eqs.gfflSl). 

Finally, we averaged over the positions of the nuclear 
spins to remove the sample dependent information. This 
revealed that the total HF Hamiltonian is quite generally 
anisotropic and, moreover, that the contribution due to 
P-like states can dominate over the contact HF contri- 
bution, see table HI1 and Eq. (l52"1) . Therefore, it can be im- 
portant to include the HF interactions (|lbllcp relevant 
for P-like states for the HESs. 

In passing, we remark that the nuclear spins can open a 
very small energy gap in the HES spectrum. This can be 
shown by averaging out all spacial directions of the nu- 
clear spin positions in the total HF interaction. Treating 
the nuclear spins as a semi-classical field of zero mean 
value ; 98 ' 99 the energy gap becomes proportional to the 
in-plane field. The ensemble averaged energy gapiSS is 
proportional to N~ x / 2 and estimated to be on the order 
of 10~ 4 /zeV for a micron sized edge, which seems out of 
the current experimental range. 
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Appendix A: On the BHZ model states 

This appendix describes various details of the BHZ 
states |P±) and \H±). In particular, the time-reversal 
properties and the phase conventions of the envelope 
functions are discussed. 

The states in the BHZ model as presented in Ref.|25| 
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are given by 

\E+) =/ 1 (z)|r 6 ,+l/2) + / 4 (z)|r 8 ,+l/2), (Ala) 

|ff+> = / 8 (*)|r 8 ,+3/2>, (Alb) 

|S-)=/ 2 (z)|r 6 ,-l/2) + / 5 (z)|r 8 ,-l/2), (Ale) 

|if-)=/ 6 (z)|r 8l -3/2), (Aid) 

similar to Eq.([5]), but without specifying any phase con- 
ventions for the envelope functions f n (z). The lattice 
periodic functions can be given as 



C/2,/5)]. Similarly 



25.72 



|r 6 , 
|r 6 , 



-1/2) 
-1/2) 



|5>lt), 

\S)\l), 



(A2a) 
(A2b) 



and 



|r 8) 3/2)=+-^[|P s >+*|P tf )J|t> > (A3a) 

|r 8 ,l/2) = -y||P,)|t)+^[|Px>+*|P w )]U), (A3b) 

|r 8 ,-i/2) = -y||p,)|;)--L[|p«)-»|p I/ 



|r.,-3/2>=-^ 



\P,) - i\Pv 



II), 



It), (A3c) 
(A3d) 



where the Bloch amplitudes \S), \P X ), \Py) and \P Z ) 
transform the same way as the well-known orbitals with 
the same names / 1 ' 101 The orbitals are connected to the 
spherical harmonics ! 69 ' 93 Thus, |r 8 ,mj) correspond to 



\j = 3/2, m-jj = 1, 



1/2) in the angular momen- 



tum representation using the total angular momentum 
J = L + S as a good quantum number, where S is the 
electron spin in the basis {t,i}- Likewise, \Te,m,j) sim- 
ply corresponds to the I = state. Note that the split-off 
band I\ with j = 1/2 and I = 1 is neglected in the BHZ 
model. Here \S) is chosen to be purely imaginary-^ and 
\Px); \Py) and \P Z ) to be real.— Furthermore, we follow 
the convention by Bernevig et al~± and Novik et aZjifi 3 - 
by using an overall opposite sigr*i°^ for the T 8 states 
in terms of the P states in Eq. (|A3[) compared to other 
authors. 69 ' 105 This sign change is not important for the 
purposes of this paper. 



W+f n {z) = E k=0 f n (z), for n = 3, 6 only. (A5) 
Here we have introduced the real operators 



T = E c (z) + —k z {2F(z) + l)k z , (A6a) 

Zm e 

W± = E v (z) + A-k z (2 l2 (z) T li(z))k z , (A6b) 

where k z = —id Zl E c i v \ is the conduction (valence) band 
edge, m e the bare electron mass, 7^2 are the Luttingcr 
parameterspS 6 - and F(z) is a real function including the 
remote T§ bands pcrturbatively^ The parameters 71,2, 
F and E c>v are different in the HgTe and CdTe layers 
of the heterostructure, which leads to the z-dependence. 
The solution of these equations will also give the en- 
ergy for that particular solution (energy band) i?k=o at 
k = (0,0). From Eq. (|A5p it follows that we can choose 
fs(z) = fe(z), which is simply denoted as /h{z) in the 
main text. Furthermore, Eqs. (|A4[) allow us to choose 
fi(z) = h{z) and / 4 (z) = h{z), which are called /_Er 6 (z) 
and fEr 8 (z), respectively, in the main text. By compari- 
son of Eqs. (|A4[) and their complex conjugates, it follows 
that we can choose fi(z) real and f±(z) purely imaginary 
as in Ref.jz|. 

Now we turn our attention to the time-reversal proper- 
ties of the states \E±) and \H±). The time-reversal op- 
erator is defined up to an arbitrary phase factor. Here 
we use O = —ia y K 7 where K is the complex conjugation 
operator and a y a Pauli matrix in electron spin-space. 
The time-reversal operator acts differently in different 
bases (due to the complex conjugation), so one should 
stick to the same basis through out a calculationi^ 3 . The 
|rj,mj) states under the time-reversal operator follow 
from Eqs. (|A2IA3|) by using that the P-like states are real, 
the S- like state is pure imaginary and that <d\f) = +\l) 
and9||) = -|t), i.e. 

9|r 6 ,±l/2) = T|r 6 ,Tl/2), (A7a) 

e|r 8 ,±i/2) = ±|r 8 ,Ti/2), (A7b) 

9|r 8 ,±3/2) = T|r 8 ,T3/2). (A7c) 



Next we discuss the phase conventions for the envelope 
functions made in the main text. The envelope functions 
fn(z) (n = 1, . . . , 6) are found from the Luttinger-Kane 
model at k x — k y — and therefore has to fulfill the 
following differential equation a 25 ' 72 

II 



Tfn(z) 

2 



P id z f n+3 (z) = E k=0 f n (z), (A4a) 



Poid z f n (z) + W-f n+3 {z) = E k=0 f n+3 (z), (A4b) 



for n = 1,2 only [i.e. only for the two pairs {fi, fi) and 



Therefore, we can now evaluate e.g. 0|-E+) by using 
Eq. (|A7p and that fET 6 is real and Jev & is purely imagi- 
nary, which gives 0|P+) = —\E—). Hence, our conven- 
tions lead to 



©|p±) = t|Pt), 

@\H±) = T \Ht), 



(A8a) 
(A8b) 



which fulfill 2 = — 1 as expected. We remark that that 
Rothe et al£2- find opposite signs under time-reversal (i.e. 
0|P±) = ±|P=F) and e\H±) = ±\H=f)), simply because 
an opposite overall sign was chosen in the definition of 
the time-reversal operator J^ 7 - In Ref . [631] . the same signs 
as in Eq. (|A8[) are found. 
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Appendix B: Normalization of the BHZ states 

In this Appendix, the normalization of the envelope 
functions and lattice periodic functions within the enve- 
lope function approximation is discussed. To this end, 
we use </?k,ff+(r) in Eq. (|7bp as an example. The entire 
wavefunction is normalized in the usual way, i.e. 



dr\ip KH+ {r)\' 



1, 



(Bl) 



where V is the volume of the entire system. The normal- 
ization of the entire wavefunction (|B1|) leaves a freedom 
to normalize the envelope function and the lattice peri- 
odic function in the most convenient way for the problem 
at hand. Various choices are found in the literature, see 
e.g. footnote 2 in the review of Coish and Baugh |28j . 

To see how this normalization choice works in prac- 
tice, we begin by separating the left-hand side of the nor- 
malization condition (|B1[) into a product of the envelope 
function and the lattice periodic function normalization, 
respectively To this end, the normalization condition 
(|B1|) is rewriting by dividing the integral over the entire 
space into a sum of integrals over the unit cells as in 
Eq.©, i.e. 

1 = -^- f dr\f H (z)\ 2 \u rg , + s(r)\ 2 

= -rrY,f, dp\f H ( Pz + z n )\ 2 \u rs ^(n n + p)\ 2 

^TtT. [ {n} dp\fH(Z n )\ 2 \u r ^ +l (p)\ 2 



J V tz Jv ^c 



L>xL>-u 



Eim z «)i 5 



K„ 



/ dp\ u r 8 ,+%(p) 



, (B2) 



where we used in the third equality that the envelope 
function - by construction - is slowly varying on the scale 
of the unit cell, so Jh{Pz + Z n ) ~ /#(Z„), and that the 
lattice periodic functions are periodic with the lattice, 
i.e. ii rg]+ 3 (7Z.„ + p) = u Tg + s_(p) for all lattice vectors 

TZ n . The integral of \u Ts+ 3(p)\ 2 over the n th unit cell 
is the same for every unit cell and hence independent 
of n, which we indicate by VuJ — > v uc . Moreover, the 
sum over lattice points in Eq. (|B2p can be made into an 
integral (including the unit cell volume v uc ), since the 
envelope function varies slowly on the inter-atomic scale, 
i-e- J2-K n \fH(Zn)\ 2 * ±-J v dK\f H {Z)\\ Therefore, 
we arrive at 



Vg 



1 



L x Ly Jy 



dn\f H (z)f 



dp\u rR , + a(p)| 5 



, (B3) 



where the normalization of the entire wavefunction in 
Eq. dBip have been written as a product of the normal- 
ization of the envelope function and lattice periodic func- 
tion part, respectively. Thus, it is now clear that some 



freedom exist in the normalization choice. 

In this paper, we normalize the lattice periodic func- 
tion as Fischer et al.r^— i-e. 



dp\u 



r 8 ,+ 






(B4) 



using that a zincblcndc crystal like HgTc or GaAs con- 
tains two atoms per unit cell, v uc = 2v a . This normal- 
ization have the advantage that the atomic HF constants 
found in the main paper are independent of the num- 
ber of atoms in the unit cell<2£ Moreover, the envelope 
function is normalized as 

-i- / dn\f H (Z)\ 2 = 1, (B5) 

L x Liy Jy 



such that Eq.(TB3l) is fulfilled. 

The normalization procedure follows the same lines 
as above for the other BHZ basis functions, e.g. all 
lattice periodic functions are normalized to the num- 
ber of atoms in the unit cell. When the wave function 
is not a simple product of an envelope function and a 
lattice periodic function, then it should be used that 
different lattice periodic functions arc orthogonal, i.e. 



■Le d P l 



,(P)«T»,mj(P) 



2Sr 



i ,!"(, Omj ,m'- • 



Finally, it 



should be noted that for the E± states, we end up with 
a combined normalization for the two envelope functions, 
i.e. J dz[\fEr 6 (z)\ 2 + \fEr s (z)\ 2 ] =1. We assume that 
each of these two envelope functions are normalized to 
one half<2ii 



Appendix C: Details on the calculation of the atomic 
integrals of the HF interactions for P-like states 

This Appendix deals with the integrals over the atomic 
wave functions of the form 






(CI) 



which appear in the matrix elements of -Hhf,2 and i?HF,3 
in Sec. IIVCI The atomic wave functions are written as 
*ffi£( r ) = i? Hg/Tc (r)Y r8im3 (0,0), i.e. a product of a 
radial and an angular part as in the main text. The an- 
gular part of the wave functions Yr g , m {Q,^) are combi- 
nations of the usual spherical harmonics^ Yi m (0, <f>) and 
the electronic spin- 1/2 (|f) and |J.)) and inherit the sym- 
metry of the bands , 44 ' 80 i.e. 



Yr„+$(M) = 


-*?(M)lt>. 


(< 


Yr 8l+ a(M) = 


-yfn°(M)lt) 


- sj\Yhe,<t>)\i), 


Yr 8 ,-*(M) = 


-sJ\Y?{e,m- 


-^\Yr\e,m) 


Yr 8 ,-*(M) = 


-Ff 1 (0,0)11), 
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which are all cigcnfunctions of J z = L z + S z (with eigen- 
value hrrij), J 2 = (L + S) 2 (with j = 3/2), L 2 (with I = 1 
due to P-states) and S 2 (with s = 1/2). To be consis- 
tent with the BHZ model, we use the same overall sign 
as Refs. 123 a nd Il03l w hich is opposite to the one used in 
e.g. Refs .l69l and ll05l (see also Appendix |A"1 Eq. (|A3|) and 
footnote 1 104 ). However, this overall sign cancel out in 
the matrix elements between P states and therefore has 
no effect here. 

To find the integrals (jCip . the spherical approxima- 
tion Eq. (|24[) is used. This is an excellent approximation, 
since most of the weight of the integrals are close to the 
atomic core. To facilitate the calculations, the HF dipole- 
dipolc like interaction for a single nuclear spin h% Eq. (|lbp 
is rewritten as (choosing the origin at the nuclear spin, 
i.e. r„ = r - R n -» r) 



Mo 



1 



1 



"47r~ 7e7> 'r 3 (l + ^)2 



(C3) 



3z 2 



i 



■Xs+u 



: (S+J-, n + S-I+, n ) 



- S-J-,n) 



® X V (Q T I O J \ 

2 \ u x ± y,7i \ u y ± x,n) 

^(s i +si ) 

2 v y z , n ' Uz y> n J 



This is written in such a way that the integrals (like 

times 



Eq. (l24l) ) consist of a radial integral over oc 



i' 3 (i+^) 



a sum of angular integrals. The terms in the curly bracket 
become the sum of angular integrals, where the space 
dependencies are seen to form spherical tensor opera- 
tors or sums thereof. Therefore, the Wigner-Eckart the- 
orem is useful to identify the integrals that are zero, see 
e.g. Ref . [93j . As an example, the element for a Hg nu- 
clear spin between the *£> r 6 r3 / 2 ( r ) states, appearing in 

the matrix element (ifkHTlHuF^l^k'HT'), is found to be 
(after some calculations) 



drr 
o Jo 

Mo 

= 4^ 7e7 -\r3 



d9sin(e)[^ T3/2 (r)}*h^l T , 3/2 (r) 



Hg 



5 t .t' 



:TflI, 



(C4) 



using Eas. (|C2IC3p and the definition (|2"5]l. The rest of 
the integrals for h% arc found similarly. 

Finally, we note that the integrals (|C1[) involving h 3 
(|lc[) are much simpler to evaluate. The integrals in the 
spherical approximation (|24[) again separate into the ra- 
dial integral (l/r 3 ) r e times a sum of angular integrals, 
which can be found by using the rewriting L n • l n = 
L z , n I z , n + §[£+,„/_,„ + £_,„!+,„] and Eq.(jU2|). 



Appendix D: Estimation of the Atomic HF constants 

In this Appendix, we estimate the atomic HF couplings 
Eqs.dTEESl), 



A Atomic 
A S,j„ 



\ Ate 



x P,j n 



2^0 
3 

Mo 

-in 



.9eMs3j„MA r |wr 6 (Rn)| , (Dla) 

■ 1 \ J« 



9efJ.Bgj n HN(Nr s ) \a jr , 



(Dlb) 



along the same lines as Fischer et al£^ These estimates 
are given in table HI Below, we go through the ingredients 
to make these estimates. 

Within the LCOA approach Eq. fTB")) , the lattice pe- 
riodic functions within a unit cell are written as a lin- 
ear combination of the two atomic orbitals. The relative 
weight between the two orbitals is related to the ionicity 
and found to bei2LI22 

a Tc ~ Vol and a Hg ^ V(U, (D2) 

which is taken to be the same for the Tq and Tg bands4^ 



Moreover, the g-factors for the various isotopes an 



,29 



.gi99 Hg = 1-01, 02oi Hg = -0.37, (D3a) 

ffl23To = -1.47, 5 i25 Tc = -1.78. (D3b) 

These are seen to vary in sign, which is the reason for 
the sign variation of the HF couplings. 

Furthermore, to estimate the HF couplings, the atomic 
wave functions \& r m c also have to be given explic- 
itly. The angular part follow the band symmetry as 
in Eq. (|C2[) . As for the radial part, we follow Fis- 
cher et alM- and approximate it by a hydrogenic radial 
eigenfunction^ 3 - R n i(r) with an effective charge eZ e g re- 
placing the actual charge of the nucleus eZ in order to in- 
clude atomic screening effects etc., i.e. Z e g < Z. The out- 
ermost electrons in Hg (Te) have the principal quantum 

number n — 6 (n = 5) such that R v (v) — Re{5).i( r ) 

and i?p 6 g(Tc) (r) = fl 6(5)]0 (r). Clcmenti et al^.^11 have 
calculated the effective charges Z e g for various atoms and 
orbitals and found that Z e g(Te, 5s) = 12.5, Z g(Te, 5p) = 
10.8 and Z g(Hg, 6s/6p) = 11.2, which obviously is much 
smaller than the bare nuclear charges eZ = 52e for Te 
and eZ = 80e for Hg. 



Using \ur 6 (R r , 



(^r 6 ,i/2) 2 K-| 2 |*f, 6 (0)| 2 with the 



hydrogenic orbital ^ 3 r (r) = R n0 (r)Yo(9, 4>) for isotope 
j, we can now give the atomic contact HF coupling as 



.Atomic _ 2 M0^ (j n n r AT ,2, | 2 [Zcff(j,s)] 2 



ira,Qn 3 



(D4) 



where ao = 4ireoh 2 /(m e e 2 ) is the Bohr radius and eo is 
permittivity of free space. 

The hydrogenic like atomic orbitals also makes it easy 
to calculate (l/r 3 ) 3 r in Eq. ([2~5l) numerically, which shows 
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that neither the nuclear length scale r c nor r max make a 
difference in practice. Hence, we can use 



1 



drr 2 \R n ^ = i(r 



1 



[Z cS (j, P )f 
3a 3 n 3 



(D5) 



to find the HF coupling A^° mic for the P states 



Therefore, now we only need one more ingredient to be 
able to estimate the HF couplings, namely the normal- 
ization constants JVr^m of the lattice periodic functions 
in the LCOA approach Eq. ([T3"l) . The normalization con- 
dition (lB4l) lead to 



rfr 



«Tc*?°, 



+ ^)±" Hg *rf 



r ~2 



N*. 



where + (— ) corresponds to Tg (Tq). First of all, we note 
that 



Nr 



Nr 



(D6) 



due to the similar form of the atomic wave functions for 
±nij, see e.g. Eq. ([C2[) . Using the hydrogenic eigen- 
states, we can therefore now numerically find the nor- 
malization constants Nr^m- Numerically, these do de- 
pendent weakly on how the Wigner-Seitz unit cell of the 
zincblcndc crystal is approximated, in contrast to the 
unit cell integrals involving hf for i = 2, 3 in Sec. IIVCI 
We have tested various spherical and cubic approxima- 
tions to the primitive Wigner-Seitz unit cell all with 
the same volume as the Wigner-Seitz unit cell, namely 
v uc = 16|d| 3 /(3v / 3), where |d| = 0.279nm is the dis- 
tance between the Hg and Te atoms in unit cell, see 
e.g. p. 58 in Ref . [3l| . Such a weak dependence is also 
found in the estimate for GaAs by Fischer et al^ From 
our various approximate unit cell calculation, we found 
that a good estimate for the normalization constants are 

(A^i/2) 2 ~ (JVr 8 ,3/2) 2 - 3.6 and (iVr 6 ,i/2) 2 - 2.7. 
Therefore, we can use approximately equal normalization 
constants for Ur 8 ,±i/2 an( i M r 8 ,±3/2j which allows for the 
introduction of a common atomic P-likc HF constant in 
Eq. ([26|) . Therefore, we now have all the ingredients to 
make the estimates with the results seen in table |U 



as 



h 

Sic h' t — — 


C k y u C k' y d + c k y d C k' y u 


(Ela) 


h 

S k y k' y ,y — *^ 


C kyd C k' y u ~ C k y u C k'yd 


(Elb) 


h 

kyky^z ~ 


C kyU C k' y u ~ c kyd C k' y d 


(Elc) 



H\J 



J 



together with the operator % yk , y = 5 [c\ yU c k , yU + c\ yd c Kd ] 
and the raising and lowering operators Sk y k' ,± — 
Sk y k' y ,x ± isk y k' y ,y for the edge state spin. Here, for in- 
stance, Sk y k/ .+ moves a particle in the state ^jj k , into 
the state \&" k and in this sense raise the edge state spin 
(while also changing the wave vector). In the case of 
k y = k' y , the edge state spin operators (|E1[) coincide with 
the usual spin- 1/2 operators^ and 3 k k is the particle 
number operator (times h/2). 

The contact HF interaction (|3T))) in terms of the edge 
state spin operators (|E1[) becomes 



nj(y) _ 



i r n iD i (k'y—ky)y n 



\ A ky,k' y + ^-ky-k'J^^kyky, 



(*») 



+ (Ai5i-A^_w)/„«3*. 



v (x n ) 



{x n 



A*-' 1 "- 1 T ir AV^"J t c, 

~ A k y -k'/-,nSkyk^+ - A_ kyk ,J + , n S kyky - 



Using the edge state spin operators (|E1[) . the HF inter- 
action (ITU) due to the P-likc states becomes 



H 



(y) _±_ 

HF > P ~ i5h 



EE 



^{k'y-ky)y„ 



B - D 

> 

B 



x A EE 



XXn 



(,v„ 



A fc n _u, I-. n 5 k k ' ,+ + A_ fc " k , I+,n£k 



(A 



\_^nl-x,n + Vnlz,n\Skyky,z 
+ (4 X y i y - ^kt-k y ) [U, fl + &A,n] % yK I , 



where 



Appendix E: Hyperfine interactions in terms of edge 
state spin operators 



Here we reformulate the HF interactions ([59]) and (pfTj) 
for the HESs along the y-axis in order to give some more 
insights into their form. Having in mind the spin-1/2 pic- 
ture of a pair of HESs discussed in sec. IVI A[ we are lead 
to introduce the non-diagonal edge states spin operators 



\A\B ZLm ^ A P,3«)> 

^4 K B D ) A p! + 3 ( B + D ) A pI\ ■ 



(E3) 



In both HF interactions the edge state spin-flipping terms 
I±,n&k k' ,=f are seen. Moreover, in the HF interaction for 
P-like states the unusual coupling I Xtn Sk y k' ,z is found as 
discussed in the main text. Note that the terms including 
the operator 3k y k' vanish to second order in k y and k' 
in the position averaging and also in the particle-hole 
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symmetric limit, see sec. IVID1 
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